Disorder-Induced Order in Quantum XY Chains 
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We observe signatures of disorder-induced order in ID XY spin chains with an external, site- 
dependent uni-axial random field within the XY plane. We numerically investigate signatures of a 
quantum phase transition at T = 0, in particular an upsurge of the magnetization in the direction 
orthogonal to the external magnetic field, and the scaling of the block-entropy with the amplitude 
of this field. Also, we discuss possible realizations of this effect in ultra-cold atom experiments. 



PACS numbers: 

Disorder is common in a wide range of physical phe- 
nomena [1]. It typically originates from inhomogeneous 
environmental conditions and other parameters beyond 
experimental control. The time scale r of observations on 
the considered system makes it possible to divide disorder 
into two categories: Annealed disorder, that changes in 
time scales much shorter than t, and correspond to fluc- 
tuations (thermal, quantum or externally driven), and 
quenched disorder, i.e. disorder that is permanent or 
does not change significantly on the scale of r. In gen- 
eral, disorder tends to reduce material properties that 
depend on regularity, such as magnetization, mechani- 
cal strength, or electric conductivity. However, even if 
counterintuitive, sometimes the opposite is true: there 
exist systems that develop order (as measured by some 
order parameter) when in presence of disorder. Vari- 
ous instances of such counterintuitive effect occur in the 
context of annealed disorder in geometrically frustrated 
systems, such as frustrated antiferromagnets. The term 
order-by-disorder was introduced by Villain et al. [2 to 
describe the effect that occurs in the Ising model on a 2D 
square lattice with ferromagnetic interactions between all 
next neighbors, except for antiferromagnetic interactions 
along the, say, even columns. Classically, such a system 
is truly disordered at temperature T — 0, and exhibits 
a large number of degenerate ground states. Thermal 
fluctuations at arbitrarily low T induce ferrimagnetic or- 
der, in which spins at odd columns order parallelly. The 
same ordering occurs in the case of quenched disorder in- 
duced by site dilution (i.e. introduction of non-magnetic 
sites) . The phenomenon of order- by-disorder persists also 
in quantum systems, although its appearance is not as 
dramatic (for a review see [3]). 

A completely different mechanism, called disorder- 
induced order (DIO), has been recently investigated by 
some of us [4], and earlier by several other authors [5]. 
This phenomenon depends crucially on two ingredients: 
the unperturbed system must have a continuous sym- 



metry, and the disorder must break it — for instance 
by promoting a single direction in space. Let us ex- 
plain this through an example. Consider a chain of spins 
with isotropic nearest neighbor interaction in the XY 
plane, which has aU{l) rotational symmetry about the Z 
axis. As it is well known, the Mermin-Wagner-Hohenberg 
(MWH) theorem [6] states that there is no spontaneous 
magnetization at any positive temperature for the bare 
(i.e. without quenched disorder) XY model for dimen- 
sions two and smaller (in the limit of infinite systems). 
The continuous symmetry of the model is pivotal for the 
proof of the MWH theorem, since it is responsible for the 
freedom of fluctuations, or, more precisely, for the den- 
sity of states of low energy excitations in the system. Let 
us add now to the system a quenched disorder in the form 
of local external magnetic fields in the XY plane, with 
random strengths and directions. The key ingredient for 
the appearance of DIO is the relation between the dis- 
order probability distribution and the symmetries of the 
physical system. In particular, if a disorder distribution 
has the same symmetry as the system, the disorder sup- 
presses the spontaneous magnetization even more: even 
arbitrarily small disorder of this form destroys the order 
at T = in 2D, and it presumably lowers the critical 
temperature of the transition in 3D |2. However, if the 
random field is aligned only along the X axis, the dis- 
order has a symmetry different than that of the system 
- the MWH theorem no longer applies, and spontaneous 
magnetization can appear precisely because of the pres- 
ence of disorder. 

The mechanism of DIO is quite general - it can be 
applied whenever there is a continuous symmetry in the 
system that can be broken with a random field with a 
lesser symmetry. Indeed, we have recently showed that a 
disordered random field with a symmetry different than 
that of the system may increase order in the direction 
where there is no disorder [4]. We considered a two- 
dimensional array of classical spins with a ferromagnetic 
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interaction in presence of a random external magnetic 
field in the Y direction. With this choice of disorder, the 
system has a remaining mirror symmetry along the Y 
axis (which can easily be broken by appropriate bound- 
ary conditions). We argued that, in this case, the disor- 
der will always induce a non-zero magnetization in the X 
axis - i.e. in the direction orthogonal to the random field. 
Similarly, the relative phase of two-component interact- 
ing Bose-Einstein Condensates can be fixed by applying 
a suitable random coupling Raman field [H]. The same 
can be done with the relative phase between the super- 
fluid order parameter of an ultra-cold Fermi gas and the 
wave function of a molecular BEC of the same species 
[9j. In this case the disorder is provided by a spatially 
random photo-associative-dissociative coupling between 
the BEC and the Fermi gas atoms. 

In this paper we report numerical studies of one- 
dimensional quantum spin chains that present disorder- 
induced order. In particular, we consider isotropic XY 
spin chains subject to a random magnetic field along one 
axis on the XY plane. We show evidence of the formation 
of spontaneous magnetization, and a transition from that 
phase to a phase without long-range order for a critical 
strength of the disordered field. Apart from the magne- 
tization, we study the behavior of entanglement. Using 
numerical methods, we are able to identify the quantum 
critical point for a particular case of staggered magnetic 
field and the critical exponents in that situation. The 
transition that we discuss here is an analogue of the one 
that can be obtained with staggered magnetic fields in 
the gapless phases of the XX Z model [10]. The main 
difference is that we mostly consider magnetic fields that 
oscillate on a much larger spatial scale, or are fully ran- 
dom. 

We envision an implementation of our ideas using 
ultra-cold atom experiments. Technological develop- 
ments have made ultra-cold gases a powerful tool to study 
condensed matter systems with an unprecedented degree 
of control - allowing the manipulation of geometry, di- 
mensionality, and even the interaction strength between 
particles (by means of Feshbach resonances) [TTI - fT5] . An 
example of the power of ultra-cold atom systems as a 
test bed for condensed matter physics was the recent ob- 
servation of matter- wave Anderson Localization [TU [TS] 
- predicted already in 1958 in the context of the metal- 
insulator transition [IB] , and for ultra-cold atoms in 2003 
|17j . These experiments, and others, illustrated an inter- 
esting difference between cold gases and usual condensed 
matter systems: because the potentials are exceedingly 
clean, and thermal fluctuations negligible, disorder in 
cold atom experiments does not occur naturally and has 
to be introduced by design. This, in turn, translates into 
complete control over the disorder distribution: e.g. by 
changing light intensity the disorder strength can be var- 
ied continuously. Typically one can use laser speckles to 
induce static unstructured randomness (pioneered by the 



late G. Grynberg [THl HH]), or superposed lasers of in- 
commensurate wavelengths to create quasi-randomness 
(for recent reviews see [20]). Moreover, one can con- 
trol nonlinear interaction strengths using Feschbach res- 
onances, and reach in this way various regimes |21j of 
the Anderson-Bose glass phase [55] . This manipulability 
is essential for deep exploration of disorder-induced or- 
der phenomena, and for observation of the phenomena 
discussed in this paper. 

MODEL DESCRIPTION 

We consider a ferromagnetic spin chain with N spins 
1 /2 in a random external magnetic field, described by the 
following hamiltonian: 

N-l N 

^ = - E i^'^'^^' + ^l^l^') - E ^^n, (1) 

1=1 4=1 

where ct^ are the a — x, y-Pauli spin matrices at site 
j, and /i* is the random field at site i. The field points 
along an arbitrary direction n inside the XY plane and 
an — n- a. Also within the XY plane, we will distinguish 
observables (and measurements) aligned about the axis 
n with a || subscript, and observables perpendicular to 
n with a _L subscript. As mentioned before, the second 
term of Eq. ([I]) does not have the same symmetry as the 
first term (which is invariant with respect to rotations 
along the Z axis). 

The relevant order parameters are the mean expec- 
tation values of the magnetization along the parallel 
and orthogonal directions: m|| = "^J- = 

(^(T^)/A^. Typically, we also consider the local mag- 
netization m|| = (cTy) and m± = {(t\), indicating which 
regions of the chain are being discussed. Both m|| and 
m± vanish as the amplitude of the external fields ap- 
proaches zero. For large field intensity, my follows the 
local direction of the field. In this case, the average m± 
is essentially zero. 

Entanglement is known to be a good predictor of quan- 
tum phase transitions [23l - l25| . Although there is a vari- 
ety of possibilities, the observation of a singularity in an 
entanglement measure most certainly implies a second 
order quantum phase transition. In order to measure en- 
tanglement we will use the block entropy S{p), defined as 
the Von Neumann entropy of the reduced density matrix 
obtained by tracing out the degrees of freedom of A'^ — p 
spins of the chain. By means of a Schmidt decomposi- 
tion, any pure state of the system can be expressed 
as 

IV') = E^J''l'/'i""')^l^^""')' (2) 

i 

where ^')} and '^')} are orthonormal 

states in the Hilbert space of the first p and last 
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N — p spins respectively. Because of the orthonor- 
mality property, in this basis it is easy to write down 
the reduced density matrix for the first p spins, pp — 
^ AilV'f '"'''l- The positive numbers are the so 

i 

called Schmidt coefficients, and give the block entropy 
Sip)^-J2\^og,\. (3) 

i 

The value of S{p) depends on both classical and quan- 
tum correlations (such as entanglement) between the two 
blocks [1, . . . ,p] and [{p+1), . . . ,N] of the TV-spin chain. 
The so called area law says that the block entropy of a 
ground state generally scales with the size of the bound- 
ary (area) of the system [25] — except at criticality, 
where there are typically logarithmic corrections. In one 
dimensional systems, the boundary of a block is constant. 
Thus, away from the critical point, entropy saturates be- 
yond a certain block size po- S{p) = S{po) for all p > po- 
Let us analyze the behavior of the block entropy of 
our system for some limit cases. For small fields, the 
XY term in Eq. ([T]) dominates, and the system exhibits 
long range entanglement. This leads to a large number of 
non-zero Schmidt coefficients — and therefore large block 
entropy. In contrast, for large amplitudes of the field, 
the second sum of Eq. ^ dominates the behavior of the 
system: the ground state is a product state with only one 
non-zero Schmidt coefficient (that must be equal to one 
because of normalization), which gives S{p) = 0. 



NUMERICAL RESULTS 
Methods and Materials 

To obtain the ground state of finite chains for arbitrary 
configurations of disorder we employ the Time Evolv- 
ing Block Decimation algorithm with an imaginary time 
evolution [^BH28j . The algorithm is based on calculating 
Schmidt decompositions at all links of the spin chain, 
which leads to describing the quantum state through a 
product of matrices. The rank of these matrices reflects 
the number of Schmidt coefficients that are retained for 
the simulations. Therefore, slightly entangled systems 
(in terms of the number of non- vanishing Schmidt coeffi- 
cients) are described accurately by small matrices, which 
leads to a large computational speedup. Strongly entan- 
gled systems, in contrast, require very large matrices to 
be described accurately. Excessive truncation of the ma- 
trices induces a breakdown of the algorithm, although in 
general one can monitor the accuracy before this happens 
— for example, by measuring the value of the smallest re- 
tained Schmidt coefficients. We also perform additional 
tests to ensure that the numerical solution does not de- 
pend on the maximum number of Schmidt coefficients. 



The TEBD algorithm used for our simulations is par- 
ticularly efficient for one-dimensional systems with on- 
site and nearest-neighbors interactions only, as is the case 
for the XY system in presence of the external random 
field. We implemented the algorithm for finite systems 
with open boundary conditions, and for infinite systems 
with a periodic Hamiltonian by imposing the periodic- 
ity of the solution. Due to the numerical complexity of 
the algorithms, we used a wide range of resources from 
desktop computers and local clusters to the Zaragoza su- 
percomputer with up to 50 parallel processors. 

The disorder- induced order effect is symmetric with 
respect to the orthogonal direction of the disordered field. 
Because of this, the original proposal used convenient 
boundary conditions in order to lift this symmetry [3]. 
For our simulations, usually it turned out to be enough 
to impose non-symmetric initial conditions for computing 
the imaginary time evolution towards the ground-state. 

Staggered field 

We begin by reviewing the case of an staggered mag- 
netic field, ft-* — (— l)'/io. Although it is not random, its 
non-uniformity will help us gain a good intuition for the 
random case. 

We observe two distinct regimes as a function of the 
magnetic field intensity (see Fig.[T]4_). For small fields, 
a finite spontaneous magnetization arises in the direction 
orthogonal to the field. On the other extreme, for large 
ho, the magnetization in this direction is zero. Inter- 
estingly, we observe that the transition between the two 
regimes is sharp, indicating the presence of a second or- 
der quantum phase transition. Our numerical estimate of 
the critical point is he = 2.915 ±0.001, which is in agree- 
ment with previous studies . As the field intensity ap- 
proaches he from below, the spontaneous magnetization 
decays according to a power law, m±{ho) ^ {l~ho/hc)^ ■ 
Our numerical analysis gives f3 = 0.125 ± 0.002. In 
Fig. [Tj3, we see evidence that at the critical field inten- 
sity of the staggered magnetization along the direction 
of the field does, indeed, show a singularity in the first 
derivative. 

In Fig. [ip we show the block entropy 6*00 for a semi- 
infinite block as a function of intensity of the staggered 
field. Near the critical point, the entropy of a semi- 
infinite block diverges as ^ : 

5oo = ^|log2(0+a. (4) 

where c is the central charge of the underlying conformal 
field theory. The factor ^ in Q appears because we 
measure entropy between two semi-infinite parts of the 
chain with only one boundary between them. Through a 
best fit to the data shown in Fig. [Tp, we obtain a value 
of c = 0.53 ± 0.05 for the central charge. 
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Figure 1: Panel A shows the ground state magnetization in 
the direction orthogonal to the staggered field depending on 
the strength of the staggered field. In B, we see the stag- 
gered magnetization in the direction parallel to the external 
field. Panel C represents the entropy of entanglement Soo for 
a semi-infinite block. Results in panel D show the entropy of 
entanglement 5*00 as a function of the logarithm of the corre- 
lation length near the critical point he — 2.915. 



For small values of the field, the entropy diverges as 
we approach the isotropic XY critical point. For larger 
values of the field intensity, entropy decays to zero, which 
is expected as the ground state becomes a product state. 
As a curiosity, for field intensities smaller than the crit- 
ical, there is a special value ho = 2^/2 of the field for 
which the block entropy is exactly zero, and the ground 
state is thus a Neel product state 
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Figure 2: Orthogonal ground state magnetization in presence 
of a regularly oscillating field. For fc = for example, we 
see the appearance of a regime with orthogonal magnetization 
for disorder amplitudes of around 1.0 and disappearance of or- 
thogonal magnetization around disorder amplitude of about 
1.5. This confirms that uni-axially oscillating magnetic fields 
can induce magnetization orthogonal to the oscillating direc- 
tion. 
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Oscillating fields 



Figure 3: Parallel magnetization in the presence of a regularly 
oscillating field. The plot shows that the spin chain magne- 
tizes according to the external magnetic field if this uni-axially 
oscillating field has a large enough amplitude. This confirms 
that the physical intuition trivially works for strong magnetic 
fields. 



Next, we focus on the case of a smooth periodic field 
such that at site i the field is /ly = hsm{ki), where 
/c ^ 1 is the wave number of the periodic field. The sys- 
tem exhibits spontaneous perpendicular magnetization 
for small, non-zero values of /i, whereas for large inten- 
sities h the parallel magnetization follows the oscillating 
field. 

Figure [2] shows the orthogonal magnetization of the in- 
dividual spins for different amplitudes of the external os- 
cillating field. Similar to the staggered field (Fig[Tj3), we 
observe two regimes of orthogonal magnetization: pres- 
ence of orthogonal magnetization for small amplitudes 
up to a given value depending on fc, and disappearance 
thereof for larger amplitudes. The sinusoidal nature of 
the uni-axial external magnetic field is translated to a 
slight variation in the strength of the orthogonal magne- 
tization. 



In Fig. [3j we see that the magnetization of the spins is 
dominated by external field when the latter has large am- 
plitudes, confirming physical intuition. The amplitude 
region of the presumed phase transition shows a strictly 
monotone increase in the tendency of individual spins to 
align with the external magnetic field. 

We show the block entropy as a function of the site 
and disorder amplitude in Fig. |4j As expected, at very 
low fields — near the isotropic XY critical point — the 
entropy grows slowly with system size, while it stabi- 
lizes rapidly for larger h. Moreover, the saturation value 
decreases with field intensity. The entropy has a local 
maximum at a field value that coincides with an abrupt 
decrease in perpendicular magnetization. 

Figure [5] summarizes our studies of the oscillating fields 
with different amplitudes. Comparing the magnetiza- 
tion in Fig. [sJ'V to the transverse magnetization shown 



5 




site nunber 



disorder anpiitude 



3.# 



Figure 4: Ground state block entropy of a partition in the 
presence of a regularly oscillating field. This configuration 
shows a minimum for amplitudes around 1.0 corresponding 
to orthogonal magnetization and a maximum for amplitudes 
of around 1.5 corresponding to the abrupt disappearance of 
the orthogonal magnetization, further indicating the presence 
of a quantum phase transition. This plot also indicates that 
the boundary effects become negligible beyond 3-5 sites from 
the edge of the spin chain. 



in Fig. [5p, we see that spontaneous symmetry breaking 
appears near the zeros of /ly . For strong enough intensi- 
ties h, this leads to the creation of a set of "islands" of 
perpendicularly magnetized spins in a sea of transverse 
magnetization. When h is weak enough, however, the 
island size R becomes greater than the distance between 
the islands ir/k, the isolated islands merge, and there is 
non-zero mj_ of definite sign everywhere (Fig. ^p). The 
periodic transition from one phase to the other can be un- 
derstood in terms of quantum phase transitions in space 
[5U| : away from the critical points, the system follows 
the local value of the field adiabatically (Fig. [5|3) and 
remains on the corresponding phase. In our system, this 
corresponds to the regions where the field is very large 
and therefore the local magnetization is in the symmet- 
ric phase = (see [3T]). However, when the am- 
plitude of the field /ly approaches its critical value, the 
local correlation length of the system, ^ ~ l/iyl"", can 
become much larger than the rate of change of the field. 

In this regime the system cannot heal fast 
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enough (compared to the change in the field), and it be- 
gins to transition from one phase to the other, forming 
an island of broken symmetry with a random sign of to_l . 
We can estimate the size R of an island by linearizing the 
magnetic field near its zero aX Iq: j/iy | « (i — io)|. At 
the boundary of the island we have the condition ^ ~ 
which writes as \hk R\~'^ ~ \R\. Thus, the size of an 
isolated island of perpendicular magnetization results 



R ~ {hky 



(5) 



This line of reasoning, based on the Kibble-Zurek mech- 
anism, allows us to also estimate how the amplitude of 



perpendicular magnetization in the islands depends on 
h and k. In a first approximation, the total magnetiza- 
tion is constant during the transition in space (Fig.[5p). 
When the local density approximation starts to break, 
i.e. at iQ ± R, the system goes from an adiabatic to an 
impulse region, and the order parameter at this point 
must "freeze". Therefore, we can estimate the ampli- 
tude of perpendicular magnetization in the island with 
the value of the parallel magnetization at the freezing 



point, 



Using that near the critical point 



|/i'fl|i/'^ [31^, we obtain 



{hkR) 



1/5 



[hky 



(6) 



We simulated an infinite system with k — 211 /N and 
N=512. We compare with the critical exponents of a 
pure XY chain, 1/5 :^ 0.14 and v ~ 0.57 [H]. The 
results, shown in Fig. [sj give R ~ with the ex- 



ponent close to the predicted value 



v+l 



0.363. For the 



amplitude of the magnetization of the islands we obtain 
again in good agreement with our predic- 
0.089. 
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Figure 5: Results for infinite system simulated with A'' = 512 
sites and periodic boundary conditions. Panel A shows 
the magnetic field /ly . In B and C respectively magneti- 
zation along the field my and spontaneous magnetization 
mx. In panel D we show the total magnetization |m'| = 
(mY)2 + (mf,)2. 



Randomly oscillating fields 

We study randomly oscillating fields produced from a 
normal distribution of mean zero and varying standard 
deviation. A large standard deviation is equivalent to a 
large amplitude of an oscillating field. Producing pseudo- 
random numbers following a normal distribution is sim- 
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Figure 6: Orthogonal magnetization of a particular realiza- 
tion of the random external field. We see that magnetization 
in the orthogonal direction of the random field is non-zero up 
to a threshold. The sudden drop to zero at a certain am- 
plitude (0.3 in this realization) indicates the phase-transition 
between the disorder-induced order state and the state pri- 
marily following the external magnetic field. 




Figure 7: Parallel magnetization for a given random field. 
Note that no abrupt change occurs in the amplitude region 
of the phase-transition (here 0.3). At unit amplitude, the 
parallel magnetization is no longer exactly following the quick 
variations of the randomly oscillating external field. 



ilar to the possible experimental realization of a disor- 
dered field using laser speckles (which has already been 
demonstrated) . Figure [6] shows the formation of islands 
of magnetization, coinciding with positive, negative, or 
alternating regions of the parallel magnetization. Since 
the external random field oscillates rapidly, the parallel 
magnetization, shown Fig. [7] cannot no longer follow the 
field exactly, even at unit amplitude. 

In close analogy to the staggered and sinusoidally os- 
cillating field, Fig. [H] shows a maximum of the block en- 
tropy at the disorder amplitude corresponding to the dis- 
appearance of the orthogonal magnetization. In fact, the 
bulk of the material no longer shows a monotone decrease 
of the block entropy for increasing amplitudes up to the 
disappearance of the orthogonal magnetization. There 
now appears a more complex structure dumping into a 
marked minimum right before the maximum at which 
the amplitude of the orthogonal magnetization vanishes. 
The contrast between these two final extremal points re- 
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Figure 8: Block-entropy for a particular realization of the 
pseudorandom external field. Note the marked minimum 
shortly before, and the maximum at the disorder-amplitude 
for which the orthogonal magnetization disappears (0.3 in this 
realization of the disorder). This apparent discontinuity in 
the block entropy also indicates the presence of a quantum 
phase transition between an orthogonally magnetized and a 
not magnetized state. 

sembles the discontinuity observed in the staggered field 
and appears to support the claim of a phase-transition 
at the corresponding disorder amplitudes. 

The fact that even for uncorrelated magnetic fields 
there appears to be (in most cases), a region of orthogo- 
nal magnetization illustrates the robustness of the effect. 
Our studies show that the disorder amplitude for which 
the orthogonal magnetization disappears now strongly 
depends on the individual uni-axial random field config- 
uration. 

Figure [9] shows the mean value over 10 disordered real- 
izations of the perpendicular magnetization. The graph 
shows a clear average presence of orthogonal magnetiza- 
tion for small amplitudes. The mean orthogonal magne- 
tization appears to be strongest for a disorder amplitude 
of approximately 30% of the XY spin-spin correlation. 

As expected, the randomly oscillating field presents a 
number of islands of different sizes because there are re- 
gions of predominantly positive, negative or oscillating 
random values. When averaging over several realiza- 
tions, we still see a clear overall induced constant order 
for small amplitudes of the external magnetic uni-axial 
random field. 



EXPERIMENTAL REALIZATION 

The proposal below combines Raman coupling to real- 
ize the spin dependent lattice with radio-frequency tran- 
sitions to induce the desired structure of the hopping 
matrices. The main idea stems from L. Mazza et al. 
[32] , and consists of creating a lattice that traps certain 
bosonic alkali atoms in all states from the lower hyper- 
fine manifold, and atoms in a certain state of the up- 
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Figure 9: Average orthogonal magnetization in the bulk of 
the spin chain. Our results show the existence of a range of 
disorder amplitudes, for which disorder-induced order occurs. 
The maximum average orthogonal magnetization is obtained 
for disorder strengths of approximately 0.3; the average is 
taken over 10 realizations of pseudorandom external fields. 

per hyperfine manifold in between the sites of the square 
lattice (the original goal of this approach was to design 
non-Abelian gauge fields in 3D lattice, details will be 
published elsewhere). 

The main ingredient is to use the magic wavelength, 
experimentally realized in systems with ^^Rb. This is 
in principle feasible with all the alkaline atoms, since 
they share the same fine structure. We propose to use 
the wavelength A such that the contribution for trap- 
ping the state |S'i/2;m-s = +1/2) coming from the states 
\P3/2;'ms = -1/2) and \Pi/2;ms = -1/2) exactly can- 
cel. In this way, the state |>S'i/2;™s' = 1/2) feels the 
potential coming from the cr_|_ polarized light whereas 
the state \Si/2','ms — —1/2) that coming from the cr_ 
polarized light (the quantization axis of ms coincides 
with the propagation direction of the circularly polar- 
ized light). One should avoid working with light atoms 
such as ^Li due to the small fine splitting. We can, how- 
ever, take heavier atoms (like ^^K, ^^K or even better 
^^Rb or ^^Rb) and eventually pump the atoms to the 
extremal Zeeman levels, that will then serve as spinless 
atoms. This approach will limit the lifetime to be < Is, 
but should suffice to observe at least some of the physics 
of DIO. For experimental realizations along these lines, 
see for example [33] . 

We start by confining the atoms in ID, say along the X 
axis. We take two laser pulses propagating in the X di- 
rection with circularly polarized light with respect to the 
X direction. Denoting by / the nuclear spin, the general 
result of Ref. [32] is that for appropriately designed laser 
fields all F = J — 1/2 states are trapped in the ID lattice 
sites, with X now the natural quantization axis. The 
optical potential for the F = I + 1/2 manifold has min- 
ima in the same lattice sites, but, interestingly, develops 
also minima for the states |F = / -I- 1/2; = I + 1/2) 
in the middle of the links in the X directions. These 



states, on the one hand, have good overlaps with the 
F — / — 1/2 states in the basic ID lattice sites, and 
obviously can serve as intermediate states for the radio- 
frequency transitions between the F = / — 1/2 atoms. 
On the other hand, the trapping potential for the states 
in the F = / + 1/2 manifold can be quite weak, so that 
tunneling effectively dominates over interactions. The 
trapping potential for F = I — 1/2 atoms can be strong 
enough to put them in the Mott insulating regime. 

In the case of Rubidium / = 3/2, one should prepare a 
large condensate in the F = 2 manifold, and then pump 
some atoms to the F = 1, F^ = —1 states. In the strong 
repulsion regime (hard bosons regime), the Hamiltonian 
for F = 1 atoms reduces to that of the XY model. The 
uni-axial random field in the XY plane can be easily 
realized using Raman (optical or RF) transitions with 
fixed phases and random strengths, as in the proposal of 
Ref. 0. 

Interestingly, the same scheme can be generalized to 
2D in a square, and even 3D in a simple cubic lattice. 
In 2D for instance, we start by confining the atoms in 
2D, say in the Xyplane. We take two laser pulses prop- 
agating in the X direction with circularly polarized light 
with respect to the X direction. Similarly, we apply two 
laser pulses in the Y direction with the circular polar- 
izations corresponding to propagation axis Y. Again, 
the general result is that for appropriately designed laser 
fields all F = / — 1/2 states are trapped in the square 
lattice sites. The natural quantization axis for them is 
now a 45 degrees axis between X and Y. The optical 
potential for the F = / + 1/2 manifold has minima in 
the same lattice sites, but, interestingly, develops also 
minima for the states |F — / + 1/2; F^ — 1 + 1/2) and 
\F ^ I + l/2;Fy = I + 1/2) in the middle of the links in 
the X and Y directions, respectively. These states have 
good overlaps with the F = / — 1/2 states in the basic 
square lattice sites, and obviously can serve as interme- 
diate states for the radio-frequency transitions between 
the F = / — 1/2 atoms. The remaining ingredient of the 
proposal are the same as in ID. 

CONCLUSIONS 

After the effect of disorder-induced order has been 
shown in classical systems, we have presented numeri- 
cal evidence that this effect also exists in quantum sys- 
tems such as quantum XY chains. The effect consists in 
the appearance of magnetization in the direction orthog- 
onal to a spatially disordered external magnetic field for 
small amplitudes of this field. The key ingredients for 
justifying this result are values of various components of 
magnetization (which are directly measurable in the ex- 
perimental scheme discussed above), and non-monotone 
block entropy. The latter cannot be measured directly, 
but its properties can be inferred from the measure- 
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ments of density-density correlations using Bragg spec- 
troscopy, noise interferometry, and/or spin polarization 
spectroscopy. Finally, let us mention that recently an 
analogue of DIO in the time domain (i.e. with time de- 
pendent perturbations) has been proposed and termed 
rocking [M] , 
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